Optimizing Optical Quantum Logic Gates using Genetic Algorithms 
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We introduce the method of using an annealing genetic algorithm to the numerically complex problem of 
looking for quantum logic gates which simultaneously have highest fidelity and highest success probability. 
We first use the linear optical quantum nonlinear sign (NS) gate as an example to illustrate the efficiency of 
this method. We show that by appropriately choosing the annealing parameters, we can reach the theoretical 
maximum success probability (1/4 for NS) for each attempt. We then examine the controlled-z (CZ) gate as the 
first new problem to be solved. We show results that agree with the highest known maximum success probability 
for a CZ gate (2/27) while maintaining a fidelity of 0.9997. Since the purpose of our algorithm is to optimize a 
unitary matrix for quantum transformations, it could easily be applied to other areas of interest such as quantum 
optics and quantum sensors. 

PACS numbers: 03.67.-a, 03.67.Lx, 42.50.Dv 



Linear optics is an attractive candidate for building quan- 
tum computers in large part due to Knill, Laflamme and Mil- 
burn and their scheme for non-deterministic quantum gates 
with projective measurement. The scheme provides a way to 
build elementary quantum gates with only linear optical ele- 
ments. The trade off in this scheme is that we can only con- 
struct the gates with a certain success probability, i.e., they are 
non-deterministic. Therefore, one of the major tasks for this 
program is to figure out a general way to determine the design 
of gates which possess the highest success probability. 

Any quantum gate or linear optical quantum state generator 
(LOQSG) can be viewed as a unitary transformation which 
transfers certain input states into desired output states. The 
goal of designing a LOQSG is to find a proper unitary matrix 
whose elements can then be implemented with linear optical 
devices fid. In this paper, we introduce genetic algorithms 
with a simulated annealing to the problem of optimizing a LO- 
QSG. We first restate the problem so that it is suitable for ge- 
netic algorithms. We then briefly discuss the method of ge- 
netic algorithms with simulated annealing and a tunable con- 
trol constraint. Using this method, we first test the efficiency 
of the algorithm with the nonlinear sign (NS) gate. We choose 
the NS gate due to its maximum success probability without 
feedforward having been theoretically proven to be 1/4 
We then investigate the CZ gate and try to obtain the global 
maximum for success probability while maintaining a high 
level of fidelity. 

Any LOQSG can be represented as in Fig. Q] Suppose we 
have /V input channels. They are composed with computa- 
tional input states and ancilla channels. We want these inputs 
to be transformed to our expected output with projective mea- 
surement on the remaining ancilla ports. This process can be 
done by a linear optical device which we call a LOQSG 
This device is an /V dimensional unitary transformation. When 
a projective measurement determines a certain pattern of pho- 
tons measured in some M < N of the modes, it is considered 
successful, which leads to a preparation of the desired state 
in the remaining modes. Therefore, the device is probabilistic 
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FIG. 1 : A Prototype of Linear Optical Quantum State Generator. It 
exploits linear operations, which eventually can be represented as a 
unitary transformation, and projective measurements to convert an 
input state into a target output state. 



and it can fail in two aspects. Firstly, the measurement does 
not give out the expected pattern, which leads to the failed 
measurement. This can be improved by increasing the suc- 
cess probability of the device. Secondly, it does not provide 
the expected output state in the computational channels even 
when the measurement works perfectly. We note this kind of 
failure as fidelity less than unity. Since we can not measure the 
computational output during the computing process, we have 
to make sure that the fidelity is equal to one or numerically 
very very close to one. 

The linear optical measurement-assisted transformation 
works as follows. We start from a computational input state 
\ijj c ) of N - M modes, combined with ancilla state \if/f n ) in M 
modes so that the input state | M^,-,, > = \if£) ® \iftf n )- The op- 
tical device induces a unitary transformation U of the IT*,-,,) 
state. After that a number-resolving photocounting measure- 
ment is applied to the M ancilla modes. The latter is formally 
described by a Kraus POVM operator in ancilla modes P = 
\vacuum A ){kN-M+\,kN-M+2, The resulting transfor- 

mation of the computational state \if/f n ) is a contraction quan- 
tum map \^o Ut ) = A|^?)/||^?|| |7[], where the action of the 
linear operator A is given by the following projection 
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In the context of the LOQSG problem, operator A contains all 
the information of state transformation. 

The optical interferometer is considered formally as canon- 
ical transformation of creation operators a] — > 1/,-,-at in- 
duced by an W x W unitary matrix U. If the input state is 
given in the Fock representation as | = \ti\,n 2 , ... , «w-m) ® 
\iin-m+\, ■•■,«/v), the unitary transformation U in equation ((TJ 
is given by 



,) = U\%„) = f] 
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Transformation of Eq. © is a high-dimensional irreducible 
representation of the N x N matrix of the optical transfor- 
mation U [8]. In Fock representation, matrix elements of 
{n\U\n') are calculated as permanents of matrix U 01 . 

Now we specify main properties of the Eq. (fl3 relevant to 
numerical implementation of the optimization algorithm. In 
the computational Fock basis \n c ), the Eq. (HJ is described 
by matrix A„c >n c = {n c l \A\n c 7 ) , which is simply a submatrix of 
{n\ U\ri). Thus A has a form of a set of polynomial functions in 
variables Uij, computed using Eq. ©, so that Eq. (0 specifies 
explicit algebraic form of dependence of A on U. If the total 
number of measured photons in ancilla modes Sjlw-M+i k I s 
the same as the number of input ancilla photons Yj^n-m+i w " 
then Eq. ([D leaves the number of computational photons in- 
variant. Since the dual-rail computational basis is just a subset 
of all possible states in the computational modes, the transfor- 
mation matrix A^jf is in general a non-square matrix, map- 
ping the Hilbert space of the computation basis to a larger 
Hilbert space. For example, A,,^ for the CZ is a 10 x 4 ma- 
trix. 

We now introduce the notion of operational fidelity of a 
transformation, which in general differs from the common 
measure of fidelity for a state transformation. From a phys- 
ical point of view, the transformation A satisfies a 100% fi- 
delity criteria if it is proportional to the target transformation 
operation A T , i.e., A = aA T , where a is an arbitrary com- 
plex number. Since the target transformation is supposed to 
be a unitary gate, i.e., Aj ar AT ar = I- The operational fidelity 
condition also requires that desired transformation A satisfies 
operational unitarity condition A 7 A - SI, where S = \a\ 2 is 
the success probability of the transformation JjJ. To formulate 
an algebraic estimate for the accuracy of the transformation, 
we consider complex rays BA and aA T as elements of com- 
plex projective space. The measure of closeness of elements 
in such a projective space is given by the Fubini-Study dis- 
tance 
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where the Hermitian inner product is (A\B) = Tr{AB^)ID c , 
and D c is the dimensionality of computational space. 

In the numerical implementation of the optimization we 
used a nonsingular variable F = cos y 2 , which we will refer 
to as fidelity in the rest of the paper. 



Success probability S of the transformation depends on 
the initial state |t^ c ) . The upper bound of S is determined 
by the operator norm ||A|| Mflx = Max((i// c \PA\if/ c }/{>p c \tJ/ c }), 
and correspondingly the lower bound of S is ||A|| Mm = 
Min{{\jj c \A' A\^/ c ) I '{\l/ c \^/ c )). As a measurement of the suc- 
cess probability, we use the Hilbert-Schmidt norm ||A|| (HS:) = 
V7>(A1A)/D C . It is easy to verify that ||A|| M/ " < \\A\\ {HS) < 
\\A\\ Max . As fidelity F -> 1, ||A|| M "7||A|| Mflv ^ 1 and S be- 
comes a well defined parameter equal to ||A|| (H,S) . We will 
refer to the Hilbert-Schmidt norm ||A|| (HS) as success proba- 
bility, keeping in mind that such a definition may not corre- 
spond to a success probability of transformation of specific 
state initial state. 

As an example, consider first the nonlinear sigh (NS) gate. 
The NS gate with 2 ancilla modes is as follows 

a|010)+yS|110) + y|210) or|010> + B\l 10) - y|210> (4) 

where \a\ 2 + \B\ 2 + \~y\ 2 = 1. Because of conservation of the 
number of photons, the coefficient matrix A has a diagonal 
form (3 x 3) with entries given explicitly as functions of U, hm . 

|010>:A11(C7) = U 22 , (5) 
|110>:A22(Z7) = U 12 U 2 i + U n U 22 , (6) 
|210> : A33(£/) = U 2 n U 22 + 2U n U 12 U 21 . (7) 



Then the operator success probability is 

S(U) = -]T\Aii(U)\ 2 . 



(8) 



Since the target matrix A^s 
is trivially calculated as 



F(E7) = 



Diag{\, 1, -1), the fidelity of A 
|Aii(f/)+A 22 (t/)-A 33 (f/)| 2 



3I? i JA H (U)\ 2 



(9) 



Similarly, The CZ gate, constructed with four ancilla 
modes, can be represented as 

or|1010> ® |1010) + /3|1001) ® |1010) 
+7101 10) ® |1010) + 5|0101) ® |1010) 
— > a|1010) ® |1010) +j8|1001>® |1010) 
+y|01 10) ® |1010) - 5|0101) ® |1010). (10) 

where + \B\ 2 + \y\ 2 + \5\ 2 = 1. In this case the matrix A is a 
10x4 matrix, and the success probability has the form 



5 ^ = iZ |A ^ )|2 



(11) 



where i=l,...,4 and j=l,...,10. The corresponding fidelity 
function is 



\Au(U) + A 24 (U) + A i6 (U) ~ A 47 (U)\ 2 
4Z U \Aij(U)\ 2 



(12) 



3 



where A 13, A24, A36 and A47 are coefficients of output states 
11010), 11001), |01 10) and |0101), respectively. 
If we write the unitary matrix as 



Comparison of Different GAs 



U = U expt^ Xigi], 



(13) 



where x is an R dimensional vector which identifies the uni- 
tary matrix U in an R-D unitary matrix space with g, as basis 
iTToh - Each point of this space is a N x N unitary matrix. If 
we take Taylor expansion on the matrix exponential term, and 
truncate the polynomial with a proper error control, then we 
can represent the U as a polynomial function of x. Substitut- 
ing it into the equations for success probability and fidelity, 
as derived in the last section, we have a polynomial function 
of x. Take x as an individual or an abstract genotype in the 
language of genetics, we then have a formula ready for the 
genetic algorithm. 

We decided to explore genetic algorithms for the follow- 
ing reasons. Firstly, the problem is to optimize a multidi- 
mensional nonlinear function (eight modes for CZ jl ill ) with 
nonlinear constraints, i.e., a very large search space. In such 
cases there is no known traditional algorithm that has proven 
promising. (See Ref. [6] for other approaches.) Genetic Algo- 
rithms (GAs) are therefore an attractive candidate. Secondly, 
the GAs can handle any form of function to be optimized, 
which will allow us to build a general scheme for our prob- 
lem. In this paper, we focus on the NS and CZ gates, but the 
method can easily be adapted to any other LOQSG device. 
Thirdly, GAs are designed for searching the optimized results 
in a global space. 

The first step in setting up a genetic algorithm is to choose 
a suitable fitness function. The simplest way is to select the 
fidelity given in Eqs. (0 and ( fTSl i as the the fitness function. 
Once we get a maximum fidelity, we can substitute the corre- 
sponding x into the function to obtain the success probability, 
Eqs. ([8]) and (fTTb . This approach is simple to implement but 
we only search the fidelity without worrying about the suc- 
cess probability. Since we know the fidelity has to be close to 
one to have a reliable result, we can consider fidelity as a con- 
straint and success probability as the fitness function. Many 
computer scientists are working on constrained genetic algo- 
rithms 1 12]. One of these is called the simulated annealing 
method 11311 . Using this method, we need to reformulate the 
fitness function so that it contains the constraint. The new 
fitness function can be written as 



(jc) = a(F(x),T)S(x), 



with 



a(F(x), T) — e 



-(l-F(x))/T 



(14) 



(15) 



where F(x) is the fidelity function used as the constraint. The 
second parameter, referred to as the temperature T, is a func- 
tion of the running time of the algorithm; T tends to (or very 
small values numerically) as execution proceeds. S(x) is the 
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FIG. 2: The success probability for NS gate of different genetic 
algorithms. It shows the efficiency for different approaches: The 
black-triangle line indicates the result of a GA without any con- 
straint, which takes fidelity as a fitness function, then calculates the 
success probability directly. The pink-square line describes the GA 
with a static-penalty when considering success probability as fit- 
ness function and fidelity as constraint, which is equivalent to set- 
ting the temperature in Eq. l !14t to be a small constant (e.g. 1CT 5 ). 
The blue-diamond line is the one which sets the annealing rate as 
T(f) = —tan~ l (t) + n/2. The green-dot line sets the annealing rate as 
T(t) = l/s/i, where t is time. 

Success Probability for CZ gate using GA 
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FIG. 3: The success probability for the CZ gate. Our goal is to ver- 
ify and exceed if possible the known maximum success probability 
2/27 ~ 0.074 which is indicted by pink-dot line. The best result 
obtained so far is 0.0737 as indicted by blue-diamond line. 



success probability, and a{F(x), T) acts like a penalty so that 
the constraint can finally be satisfied. When the GA begins, 
we want the penalty to be small, i.e., a » 1, so that the algo- 
rithm can search a bigger space to find the global maximum. 
When T is large, which happens at the beginning of the exe- 
cution, then at « 1. As time goes on, T — > 0, then a — > 0. 
It means the fitness tends to zero unless the constraint is sat- 
isfied, i.e., F(x) m 1. Therefore, at the end of the GA run, 
we can get the optimized success probability with a fidelity of 
one. The details of this simulated annealing genetic algorithm 
are described in Ref. [13]. 

The annealing genetic algorithm provides a way to search 
the global maximum of success probability of a LOQSG sys- 
tem and guarantees that fidelity is equal to one at the same 
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time. The disadvantage of this approach is that the efficiency 
delicately depends on the choice of the temperature annealing 
rate 111 211 . In the following section, we investigate this problem 
using the NS gate as an example. 

We use EOlib [ 14] as the genetic algorithms framework. It 
provides a basic genetic operation template. In the case of 
the NS gate, we compared efficiency for different approaches. 
The comparison is shown in Fig. [2] The vertical axis de- 
notes the success probability and the horizontal axis repre- 
sents the events corresponding to each implementation. Each 
point indicates a complete run starting from a randomly se- 
lected population (a set of x. We rearrange the points in 
ascending order so that we can compare the efficiency eas- 
ily. In Fig. [2] The black-triangle line indicates the result of 
a GA without any constraint, which takes fidelity as a fitness 
function, then calculates the success probability directly. The 
pink-square line describes the GA with a static -penalty when 
considering success probability as fitness function and fidelity 
as constraint, which is equivalent to setting the temperature 
in Eq. (fT4l to be a small constant (e.g. 1CT 5 ). The third line 
with blue diamonds is the one which sets the annealing rate as 
T(t) = —tan~ x (f) + n/2. The fourth green line sets the anneal- 
ing rate as T(f) — 1 / y[i, where t is time. From these plots, 
we can see that there are more chances for genetic algorithms 
with an annealing penalty to reach the global maximum of 1/4. 
The most impressive result is that when choosing the anneal- 
ing rate as 1 / yfi, the simulation can get the global maximum 
for each run. This is also consistent with the practical estima- 
tion inRef. fll . 

Using this approach, we now design a CZ gate with high 
success probability. The results are shown in Fig. [3] They 
show a strong support for Knill's highest known success prob- 
ability (2/27) 1 1511 being the actual maximum success proba- 
bility. Due to the complexity of the CZ gate, it is not known 
if an analytical proof for determining the maximum success 
probability is possible. As discussed, once we get the x, we 
can construct the corresponding unitary matrix using Eq. ( TT~3T >. 
We can then discover the optimized design for the quantum 
circuit. 

In this work, we introduced simulated annealing genetic al- 
gorithms to look for an optimized linear optical quantum state 
generator. By investigating the NS gate, which has been the- 
oretically proven to have a maximum success probability of 
1/4, we found that we can reach the global maximum for each 
run if we carefully choose the annealing rate. We have found 



what we believe to be an upper bound for Success Probability 
for the CZ gate, our results showing a plateau around the best 
known published results of 2/27 » 0.074 |15|]. Our best result 
so far having been 0.0737. Since this approach is focused on 
searching the unitary matrix space to optimize the quantum 
circuits, it can be generalized to other devices, for instance, to 
quantum sensors. 
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